On the Quantum Kinetic Approach and the Scattering Approach to Vacuum Pair 

Production 



Cesim K. Dumlu 
Department of Physics 
University of Connecticut 
Storrs, CT 06269-3046, USA 
E-mail: dumlu@phys.uconn.edu 

For Schwinger pair production with background electric fields varying only in time, it is shown that 
the quantum kinetic approach, based on the quantum Vlasov equation, and the quantum mechanical 
scattering approach, are completely equivalent. This different perspective on the quantum kinetic 
approach suggests new approximation methods, which are tested in detail using a soluble pulse-like 
field configuration, and may give useful insight for generalization to more general electric fields, 
varying in more than 1 dimension. 



I. INTRODUCTION 



It is well known that the quantum loop corrections to the classical Maxwell Lagrangian lead to the nonlinear 
interactions of light due to vacuum polarization [l|, Q ■ Some of these effects such as the Casimir effect and nonlinear 
Compton scattering have been experimentally observed. Others such as vacuum birefringence and electron-positron 
pair production from vacuum have still to be directly observed. Definite plans exist for such observations [E HI- At 
planned laser facilities, such as the Extreme Light Infrastructure (ELI), extreme field strengths might be reached and 
the phenomenon of pair creation in an external field could be observed 

The crucial goal for theory is to find a quantitative relationship between the applied electromagnetic field 
and the rate at which vacuum pair creation occurs. The problem formally can be worked out using quantum field 
theory as the pair production rate can be deduced from the imaginary part of the effective action . However clear 
explicit results exist only when the electric field is a function of a single coordinate variable. If this is the case, 
more straighforward approaches exist to find the pair production rate. For example by representing the laser field 
as a time dependent but spatially uniform electric field, one way is to reduce the time evolution to a 1 dimensional 
scattering problem [E 0], which can then be solved using the numerical or the WKB methods [E 0, H, IE EE]- 
Another is the quantum kinetic approach in which the adiabatic particle number satisfies a quantum Vlasov equation 

In this article, we show that these two rather different-looking approaches arc in fact equivalent. This new 
perspective leads to new and improved approximation methods. 



II. METHODS AND EQUIVALENCE 

A. Quantum Kinetic Approach 

In the quantum kinetic approach, the adiabatic particle number operator for scalar and spinor QED is defined in a 
time dependent basis to realize the effect of the time dependent background field [ll|, EE Ellil EE El- The adiabatic 
particle number is then defined as the expectation number of this operator: 



McW = <4(*K(*)> = (&l k (f)&-k(t)>, (i) 

where a_k(t)(aL k (i)) and b-k(t)(b_ k (t)) correspond to the annihilation (creation) operators for the particle and 
the antiparticle respectively with the momentum |k|. The expectation value can be labeled by momentum since the 
background electric field is homogeneous in space. To specify the time dependent basis for the creation and annihilation 
operators, for scalar QED, one starts with the Klein Gordon (KG) equation. Using the Fourier decomposition of the 

field operator, §(x, t) = J e l ^ ,x ^Mk(i) a k + u kW^-k) > ano - ^ nc ^ ac ^ that the gauge field varies only in time, the Klein 
Gordon equation reduces to: 



d 2 u k (t) 
dt 2 



ul(t)u k (t) = 0. (2) 



2 



Motivated by the constant electric field case, one can make the adiabatic ansatz for u\^(t): 

Uk ( t ) = a ^ e -iMt) + gkW e i M t) (3) 

where a k (i) and /3k(*) are the Bogoliubov coefficients. Using the gauge A(t) = —E(t), the functions u> k (i) and (j»k(t) 
arc defined as: 

Wk(«) = Je\ + k\{t), M*)= I ^(t')dt', (4) 

fe|| (t) = fc|| - ?A(f), e ± = y/m? + kl, (5) 

where fen is the momentum along the direction of the applied electric field. Consistency of the ansatz (3) leads to a 
coupled set of differential equations for the Bogoliubov coefficients: 

A similar decomposition exists for the spinor QED, based on the Dirac equation rather than the KG equation. In a 
suitable spinor basis the time dependence of the field operator can be expressed again in terms of a scalar coefficient 
function Uk(t) satisfying fl6| : 



dt 2 

Making the adiabatic ansatz: 



+ (w£(t) + ifc||(t))uk(t)=0. (7) 



uu(t) = P -^ k (t) , e i^{t) ( 8 ) 

^2uj^{t){uj^{t) - k\\{t)) ^ k (t)(^kW-fc||(t)) 



a k (*) -i4> k (t) Pk(t) 



y2w k (t)(w k («)-fc||(t)) ^2w k (i)(wk(i)-fe||W) / 

leads to the coupled differential equations for the Bogoliubov coefficients: 

The transformation between the time dependent operators (ak(i), &k(i)) and the time independent operators (cik, &k) 
is expressed in terms of the Bogoliubov coefficients a-^(t) and Pk(t): 



o k (t) 



(10) 



Here the upper/lower sign refers to scalar/spinor QED, respectively. Imposing the field operator commutation rela- 
tions, the relationship between the Bogoliubov coefficients can be given for scalar (upper) and spinor (lower) QED 
respectively: 

KW| 2 T|/3k(i)| 2 = l. (11) 
Then, for both bosons and fermions, one finds: 

A/k(*) = IAcWI 2 . (12) 

The number of pairs produced with the momentum k is given by the asymptotic value A/k(oo). The quantum kinetic 
equation for Afk(i) follows from the time evolution of the Bogoliubov coefficients in (6) and (9). We can express the 
equations (6) and (9) as (from now on, the subscript k is suppressed for notational simplicity): 

a ± (t) = ^l(3 ±{t )e^\ p ± {t) = ^MatWe-^), (13) 
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where the functions W±(t) (for scalar and spinor QED) are defined as: 

W + (t) = *$, W _ {t) = ^-. (14) 

Taking the time derivative of (12), and using (13), one obtains the equation for the evolution of the adiabatic particle 
number. Further with the aid of (11), one gets the quantum kinetic equation which is also known as the quantum 
Vlasov equation [H GJ, El EI El : 

Af±(t) = W± ^ ^ J W±(t)(l ± 2A4(i'))cos (2 J drw(r)dr^ dt'. (15) 

The important physical characteristic of the quantum kinetic equations is their non-Markovian nature. The change in 
the particle number at a moment t depends on the particle number's past configuration. In other words, the integral 
incorporates memory effects. 



B. Scattering Picture 

Another common approach to the vacuum pair production problem when the electric field just varies in time, is 
based on a WKB approximation to an effective one dimensional quantum mechanical scattering problem. In this 
approach [f| 0, E3, EH , equation (2) is interpreted as a one-dimensional Schrodinger equation with t playing the role 
of x, and with "potential" given by uj 2 {t). The physical boundary conditions for vacuum pair production correspond 
to the scattering boundary conditions, and the reflection coefficient can be related to the adiabatic particle number. 
Often this reflection coefficient is computed using the WKB approximation, but it can just as well be computed 
exactly (numerically). For example for scalar QED, Popov proposed the ansatz Q: 

u{t) = B{t){e-^ {t) - R{ty^), (16) 
ii(t) = iu(t)B(t)(-e- lc t>^ _ ij(t) e <*(*)). 

The second line requires B(t) and R(t) to be related through: 

B(t)( e -**W + i?(t)e^ (t) ) + R(t)B(t)e l4,{t) = 0. (17) 

This ansatz converts the "Schrodinger" equation (2) into a Riccati-type equation for R(t) 

R( t ) = (V 2 ^» - R\t)e 2i ^) , (18) 



with initial coefficient condition R(— co) = 0. Comparing (16) and (18) with (3) and (9) respectively, we identify 

m 



R(t) = ^ftj, and |i?(oo)| 2 corresponds to the reflection coefficient. Thus, 



| M2 _M ATM 

|ii(oo)l ~ I^MF ~ l+N{oo) (19) 

A similar argument for spinor QED, based on (8) and (9) can be given. Together with scalar case the Riccati equations 
are given as: 

±W ± (t) f„_ M<Kt) R 2 (i)p 2 l4 >(t)\ , D ,_M2 _ A4(00) 



R±(t) = —f^ [e-«+w T Ri(t)e^>) , 1^(00)1" = 1± ^ ±(w) (20) 

This argument shows that the quantum kinetic theory evolution equations (15) are equivalent to the 1 dimensional 
scattering Riccati equations (20). This is not immediately obvious looking at (15) and (20). To see this equivalance 
more explicitly, we note that the Riccati equation (20) involves both the amplitude and the phase of the reflection 
amplitude R(t), but only the amplitude |i?(oo)| 2 is needed. We can eliminate the phase dependence to obtain second 
order nonlinear differential equations for |i?(£)| 2 , which are naturally expressed as: 



= (*) = ^^Q±{t)±Wl{t)Q ± (t)-2uj{t)W ± (t)^±{Qlit)-l)- (w±Q±{t))\ (21) 
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where we have defined 

^ = \rww = 1±w±it) - (22) 

The equivalence of the Riccati equation (18) to (21), for scalar QED, can be shown by defining R+(t) as: 

R+(t) = M(t)e 2tx{t \ (23) 

where both M and x are reai - Time derivative of (23), with the aid of equation (18), yields a set of coupled differential 
equations: 

Mt) - W+{t Ko Smt) + 2xm (24) 



1-M 2 (i) 2 
mt)M W+{t) 
l + AI 2 (t) ' 2 



frin(20(t) + 2x(t)). (25) 



These coupled differential equations can be shown to equivalent to (21) using the definition (22) of Q(t). To show that 
the quantum Vlasov equation (15) is also equivalent to (21), we differentiate (15) and use definition (22) to obtain: 

= ^r^Q+W + KQ+it)- 2 "® I W+(t')Q + (t')sm(2 ( j ) (t)~24(t / ))dt', (26) 

W+(t) J-oc 

We recognize the first two terms in (21) immediately. The final, square root, term in (21) by using (24) and (25) can 
be given as: 

2^)^(4)-! - (w + (t)Q + (t)) 2 = \ M(t ^ sin (2<j>(t) + 2 X (t)). (27) 
To show the complete equivalence, setting last term of (26) equal to the righthand side of (27) requires: 



d(cos2x(t)jQ 2 + (t)-l 
W+(t)Q+{t) cos (20(t)) = 1 



d(sm2x(t)^/Q 2 + (t)-l^ 



W+(t)Q+(t) sin (2c/>(t)) = * ^ '-. (29) 

Above equations can be shown to be true after differentiating the righthand sides and utilizing (24) and (25) to 
eliminate xW- A similar argument can also be applied for the spinor case. In these differential equations the memory 
effects of the Vlasov equations (15) are registered in the nonlinear square root terms in (21). 



III. APPROXIMATIONS 



For realistic potentials that might realize pair creation, the Schrodinger equations (2) and (7) cannot be solved 
analytically. This fact makes it clear that numerical methods have to be used. Also note that the quantum kinetic 
approach gives jVk(oo), the number of particles of momentum k, and the scattering approach gives |i?k(oo)| 2 , 
which determines A/k(oo). Therefore, to obtain the total particle number one has to integrate the numerical results 
over the momenta, which is computationally intensive. Hence, the efficiency of these two equivalent approaches is vital. 

Several approximations can be made to improve efficiency without giving away too much from accuracy. Two 
basic approximations that can be done to Vlasov equations are the Low Density Approximation (LDA) and the 
Markovian Approximation (MA) [l6l |. It is important to stress that these approximation methods are peculiar to 
the quantum kinetic approach. In the LDA, it is assumed that the particle density J\f±(t) is very small. Then the 
quantum kinetic equations reduce to: 

Nt DA {t) - j W±{t')cOS [2 dTu(T)dT^Jdt'. (30) 
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Another approximation is the MA, in which one neglects the memory effects induced on the particle number, so that 
the factor (1 ± 2N±(t)) can be taken out of the integral, yielding: 

Afl IA (t) = -^±^(1 ± 2N¥ A (t)) f W±(t') cos (2 f dT0j{T)d T \ dt'. (31) 



These two approximations are clearly related to one another by: 

Nl IA {t) = (1 ± 2M¥ A {t))ti± DA {t). (32) 

The LDA is physically relevant when the intensity of the applied field is very low, since the pair production rate 
for low intensity fields is small. This type of weak field (qE(t) -C e]_) also ensures the applicability of the MA. This 
becomes more evident, if one classifies the relevant timescales for the pair production process. The characteristic 
timescale for particles is expected to be order T part i C ie ~ 77; , while the relevant time for pair production can be taken 
as T pair ~ q E(t) • Then from the weak field condition above, one has: 

Tparticle < ^ ^~pair • (^^) 

The non-Markovian behaviour states that the number density configuration at time t depends on the external field's 
past configuration up to time t. The relevant timescale for particles is the measure of the time required to induce 
memory effects on the number distribution. If this time is very small compared to the relevant time scale of pair 
production, which depends on the energy of the field, then non-Markovian effects cannot be resolved by the system. 
Therefore the above inequality ensures the Markovian assumption can be applied. However, this does not mean that 
the system becomes completely Markovian, since there will still be memory effects represented by the integral in 
(15), even if the (1 ± 2A/±(£)) term is taken out. 

The equivalent quantum mechanical scattering interpretation suggests other approximations. The most obvi- 
ous one is to make a WKB approximation to the scattering problem. Since the "potential" — Lo 2 {t) is negative this 
is an "over-the-barrier" scattering problem, which is conveniently described using the "imaginary time" method. 
Defining 2;' = it we find 



d 2 u{x') 



+ u 2 {x')u{x) = 0, (34) 



dx' 2 

so we have a one dimensional quantum mechanical tunneling problem. In this case, the asymptotic value of the 
particle number corresponds to the transmission or tunneling probability. Identifying the problem as a tunneling 
pro cess enables us to make the WKB approximation. In WKB approximation, the tunneling amplitude is given by 

|i?M| 2 « e - 2S *—\ ^classical = £ C(t)dt, (35) 

where ^classical is the classical action. The WKB approximation is valid in the quasiclassical regime in which the 
following inequalities hold: 

£<Cm 2 , -<m, \^.E. (36) 

Here t stands for a characteristic width of the time dependence of the external field, and E is the peak magnitude. 
The WKB approximation requires the computation of the contour integral expression for the classical action. 

Returning to the Riccati equation (20), a simpler approximation method can be introduced. From (19) it 
can inferred that in the LDA the limit Af±(i) — > implies taking the limit R±(t) — ► 0. Then, the Riccati equation 
for scalar and spinor cases can be reduced to: 

^(^Sle-^W, (37) 

neglecting the nonlinear term. We call this the Truncated Riccati Approximation (TRA). Although this approach 
seems to be mathematically equivalent to the LDA, its numerical evaluation yields almost the same results with that 
of the MA. 



A comparison between these various approximations provides valuable insight to the pair production problem. 
This can be done for a pulse shaped time dependent potential for which the exact solution is available. Therefore, a 
comprehensive assessment can be done on the basis of accuracy and efficiency. 
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FIG. 1: Quasiclassical regime, scalar case. Four approximation schemes are indicated. The plot shows the bell-shaped depen- 
dence of pair production rate on fell, in units m. The field intesity is taken as E = 0.1E cr = O.fra 2 , and the field width is set 
to r = fOm. Mass and charge are taken as f . Normalized units (h = 1, c = f ) were used for all numerical calculations. 

IV. COMPARISON 

In this section we compare the mentioned approximation methods using an exactly soluble pulse-shaped field: 



A(t) = -Ert^ t -, E(t) = J -l w , 



(38) 



for which the KG and Dirac equations arc exactly soluble Therefore, a closed analytic form of the particle 

number, for both scalar and spinor QED; can be given respectively as [7l. [l9|: 



|Af(«>)| 

|i?-M| 2 



2 cosh (a — (3) + cos A 



cosh (a + P) + cos A ' 

sinh ( 2=§B1 ) + sinh ( S^&=>L ) 



sinh ( 



sinh 



'a+P-X' 



The parameters are defined as: 



a = TTT^/ni 2 + k 2 L + {k\\ + eE) 2 , [3 = -KT^Jm 2 + k 2 L + {k\\ - eE) 2 , 

A = 



1 



\l^-{2eE T ) 2 , 



(39) 
(40) 

(41) 
(42) 



The exact closed form of the particle number becomes useful to check for computational errors and the validity 
of the numerical methods. For both scalar and spinor cases the agreement between the Riccati, quantum kinetic 
equations and the closed form is well established up to 7-8 significant figures in Mathematica for general values of 
the parameters. Here we present the results of the comparison between the approximation schemes in the relevant 
domains of momentum and electric field intensity. 

A. Scalar QED 

The first relevant domain is the quasiclassical regime where the inequalities in (36) hold. In this regime, the 
WKB approximation becomes valid and all four schemes can be compared. For the chosen gauge field, the WKB 
approximation yields the pair production rate as [Io| : 



-2S c l ass ical 



-(0+/9-Y) 



(43) 
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FIG. 2: The high intensity field regime, scalar case. Again the bell-shaped dependence of pair production rate on ku is observed. 
Electric field intensity is taken as E = 0.8-B C r = 0.8m 2 . The field width is taken as r = 10m. Mass and charge are taken as 1. 
It is observed that the rate changes drastically as the field intensity increases. 



In Figure 1. the approximations are compared with the exact closed form (39). The numerical results show that in 
the quasiclassical regime the WKB approximation overshoots the exact result when the momentum gets close to 
zero. The same pattern is observed for all remaining approximations. However they differ from the exact result up 
to ~ %5. In this regime, the TRA, MA and the LDA results are indistinguishable. 

The other relevant regime is the domain where the electric field intensity gets close to the critical field [E = m 2 , 
see Figure 2.). In this domain, the TRA and the MA give closer results to the exact answer than the LDA. However 
results of these 3 approximations lie below the analytic result unlike, in the quasiclassical regime. 



B. Spinor QED 

In the quasiclassical regime of spinor QED, the WKB approximation gives the closest result to the analytical 
expression. The closed form of the WKB method is given as [10( | : 

e -2S cla3Blcal = e _(a+/3-A') e § (V^i+V^i- 2 ) i (44) 
a = 2-KeET. (45) 

For spinor QED, the momentum distribution falls off faster than in the scalar case. However, similar to the scalar 
case, the TRA, MA and the LDA give results on top of each other. When the electric field intensity gets higher, the 
WKB approximation still gives the best result. Like in the scalar case, in this regime the distinction between the 
TRA, MA and the LDA gets clearer. The LDA differs from the other approximations by up to ~ %5. The difference 
between the analytic results goes up to ~ %25. 



V. CONCLUSION 



For intensities lower than the critical field the number of produced pairs can be calculated within the QED effective 
action formalism. For external electric fields varying only in time, the Riccati and quantum kinetic equations can 
be used for numerical evaluation and they have been shown to be equivalent. The scattering approach to the pair 
production problem gives a new perspective and suggests other approximation schemes which have been compared 
with the Low Density Approximation (LDA) and the Markovian Approximation (MA) of the quantum kinetic 
approach. The numerical data suggests that the Truncated Riccati Approximation (TRA) gives the same results 
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FIG. 3: Quasiclassical regime for the spinor case. As in the scalar QED example in Figure 1, we take E = 0.1E cr = 0.1m 2 , 
and r = 10m. Plot shows the values for pair production rate differs slightly than the scalar case. 
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FIG. 4: High intensity field regime, spinor case. As in the scalar QED example in Figure 2, we take E — 0.8E cr , and r = 10m. 



as the (MA) of quantum kinetic approach. Among the discussed approximation schemes, the TRA appears to be 
the simplest one for numerical calculation. The particle number can also be calculated by using the second order 
nonlinear differential equations featuring the non-Markovian behaviour. However, the efficiency of these methods 
requires further investigation. 

Although these two methods are equivalent, they can have different physical interpretations. In the kinetic 
approach, the pair production rate is physically related to the particle number. On the other hand, using the Riccati 
equations, the pair production rate can be associated with the reflection coefficient. Another physical interpretation 
can be made on the basis that (2), in fact, represents a harmonic oscillator equation with a variable frequency. Then 
the problem of pair creation for the scalar particles becomes equivalent to the parametric excitation of an oscillator 
[l8| . In this respect, the excitation parameter becomes the pair production rate. As for the spinor case, the same 
procedure starts with the Dirac equation and in one dimension as suggested in Q, the problem can be viewed as 
a problem of spin precession. These arguments suggest that looking at the problem from different standpoints, 
one may calculate the rate and make approximations relevant to the particular perspective taken. However, the 
applicability of the quantum kinetic approach is restricted to the fields varying in one dimension. If the scattering 
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approach can be generalized for the fields varying in more than one dimension and for one dimensional fields with 
more complicated structure [20] , it might give useful insight to the vacuum pair production problem. 
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